Coherent electronic transfer in quantum dot systems using adiabatic passage. 
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We describe a scheme for using an all-electrical, rapid, adiabatic population transfer between two 
spatially separated dots in a triple-quantum dot system. The electron spends no time in the middle 
dot and does not change its energy during the transfer process. Although a coherent population 
transfer method, this scheme may well prove useful in incoherent electronic computation (for ex- 
ample quantum-dot cellular automata) where it may provide a coherent advantage to an otherwise 
incoherent device. It can also be thought of as a limiting case of type II quantum computing, where 
sufficient coherence exists for a single gate operation, but not for the preservation of superpositions 
after the operation. We extend our analysis to the case of many intervening dots and address the 
issue of transporting quantum information through a multi-dot system. 

PACS numbers: 73.23.Hk, 73.21. La, 73.63.Kv, 03.67.-a 



INTRODUCTION 

Practical devices based on coherent quantum mechan- 
ical properties are widely seen as being one of the next 
major advances in electronics. This search is typified by 
the keen interest in quantum computation P| . In this pa- 
per, however, we seek first to explore a slightly different 
aspect of coherent quantum devices, by investigating how 
coherence can be used to enhance the capabilities of an 
otherwise incoherent electronic device. We will then con- 
sider extensions where quantum information, rather than 
classical information, is transported around a quantum 
dot system. 

We consider how to utilize coherent charge transfer 
in an incoherent device, i.e. one where the initial and 
final states of the device are eigenstates of the system 
Hamiltonian, but intermediate states are formed by co- 
herent superposition states. Applications which might 
benefit by such coherent charge transfer could be sin- 
gle electron pumps pi] or electronic digital logic (for ex- 
ample through quantum-dot cellular automata (QDCA) 
13), but our ideas are general and should have wide ap- 
plicability. By exploring this regime we seek to increase 
our understanding of quantum coherence and exploit co- 
herent effects before the technical difficulties of building 
truly global coherent quantum devices have been solved. 
As such one may even consider our proposal as realising 
a limiting case of a type-II quantum computer Q where 
there is enough coherence to effect a single gate operation 
(a swap gate) but not enough coherence to maintain any 
superpositions after the initial transfer is completed. 

Mechanisms of coherent population transfer can be 
broadly classified into two classes, nonadiabatic and adi- 
abatic. Nonadiabatic methods are usually typified by 
control of tunnelling to realize tt or similar pulses, which 
transfer population from one state to another coherently 



and rapidly. Nonadiabatic processes require exacting 
control over pulse areas and tunnelling rates that are 
fast compared with all relevant dephasing times. Al- 
though nonadiabatic population transfer is a goal, the 
sensitivity to noise of nonadiabatic transfer can place pro- 
hibitive requirements on pulse switching times in certain 
implementations. Adiabatic transfer is usually typified 
by quasi-static population transfer. Such methods are 
usually slow and are accompanied by high fidelity and 
relative insensitivity to gate errors and other external 
noise. Between these regimes lie adiabatic passage tech- 
niques which combine the advantages of robustness from 
adiabatic techniques with well-defined trajectories to re- 
alize population transfer rates approaching those of the 
nonadiabatic regime. We differentiate adiabatic passage 
from other adiabatic transfer processes by the simulta- 
neous modulation of at least two system parameters to 
realize a desired trajectory through the Hilbert space. 
For a further discussion of adiabatic and nonadiabatic 
timescales in quantum electronic transfer see Ref. 0. 

Adiabatic passage techniqes were originally developed 
for NMR and optical applications. For a good review see 
Ref. y . There have been several investigations of the ap- 
plications of electromagnetically mediated adiabatic pas- 
sage in solid-state quantum computing, see for example 
Ref. There are reasonably few proposals for adia- 
batic passage without using electro-magnetic radiation 
in the literature, although Brandes and Vorrath consider 
adiabatic passage through a double-dot system j^; Bar- 
rett and Milburn considered using adiabatic passage as 
a method for determining dephasing rates in a qubit j^; 
and more recently Zhang et al. [lOj have proposed a 
scheme for generating entanglement in a triple-dot, two- 
electron system via adiabatic passage. Adiabatic passage 
has been applied to the problems of rea ding out a charge- 
superposition qubit 'i^ and spin qubit 'll"|. 
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FIG. 1: (a) Schematic of the triple donor system realized 
with individually placed donors in silicon. Two of the donors 
are assumed ionized, the other neutral. The donors are la- 
belled from left to right 1,2,3. The position eigenstates of 
the electron are |2), |3), corresponding to the donor la- 
bels, with energy Ea for a = 1,2,3 which can be controlled 
by the S gates. Coherent tunnelling is between nearest neigh- 
bors only, with rates fii2 and 0.23- These rates are controlled 
by external B gates and can therefore be time varying, (b) 
Another possible implementation in a 2-D electron gas. Light 
gray rectangles represent the S and B gates. The gates are set 
so that only one electron is allowed in the system, with only 
one electronic state per region. There are electron reservoirs 
to the left and right hand side. Readout for both schemes is 
assumed via sensitive electrometers (not shown). 



We propose adiabatic passage through a triple- well sys- 
tem. We term this method CTAP (Coherent Tunnelling 
by Adiabatic Passage) and it constitutes an all-electrical 
analogue to the more well-known optical STIRAP (Stim- 
ulated Raman Adiabatic Passage) scheme for population 
transfer in a three-level A system. Independantly, Eckert 
et al. |l2l | have considered the realization of an analogous 
STIRAP style scheme with a neutral atom in a triple 
well potential formed in an optical trap. Such schemes 
are more complicated than two-state adiabatic fast pas- 
sage, but have the advantage that the electron energy is 
unchanged throughout the population transfer. This con- 
sideration is important when considering the trade-off be- 
tween heat dissipation and gate speed in QDCA 0], and 
our scheme could be used to realise true adiabatic clock- 
ing of a QDCA, as opposed to the quasi-adiabatic clock- 
ing schemes reported recently |l4j . Maintaining constant 
energy is also important for preventing an extra dynam- 
ical phase appearing when quantum information is being 
transported, as will be discussed below. 



TRANSPORT IN A TRIPLE-DOT SYSTEM 

A schematic of our system appears in Fig. ^ We have 
a triple- well system |1), |2), |3), where we wish to achieve 
coherent population transfer from |1) to |3) without any 



population being in state |2). The system is controlled 
by shift gates, Si which control the energy of dot i, and 
barrier gates, which control the tunnelling rate be- 

tween dots i and i + 1. Our scheme is applicable to any 
system where the coherent tunnelling times can be var- 
ied from much slower than the Ti time (population relax- 
ation) to several times faster than the T2 time (coherence 
relaxation). An ideal system to consider for implement- 
ing this proposal would be the the phosphorus in silicon 
donor scheme introduced in [isj for quantum computing, 
however it could equally well be realized in a supercon- 
ducting system 0|, two-dimensional electron gas [l^ . 
or even conceivably in a metallic triple-dot structure re- 
alised via small buried phosphorus clusters ll8|(25||. To 
proceed, wc write down the Hamiltonian for the triple- 
well system in matrix form, with state ordering |1), |2), 
|3) 
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(1) 



where flaf3 — ^ap{t) is the coherent tunnelling rate be- 
tween position eigenstates |a) and |/3) and A = E2—E1 — 
E2 — E3. In practice, to maintain the state energies will 
require the use of S gates for tuning and compensation. 
The eigenstates of the Hamiltonian of Eq. ^ are 

11?+) = sineisine2|l) +cose2|2) +00361 sine2|3), 
P_) = sin6i COS82II) — sin62|2) + cos9i cos02|3), 
\Vo) = cosOill) -H0|2) -sinei|3), (2) 

where we have introduced 

81 = arctan(f7i2/il23) , 
1 



~ — arctan 
2 



(v/(2fir!i2)2 + (2ftr!23)2) /a] . (3) 
The energies of these states are: 



£0 = 0. 



A 1 



- ± -v/(2ftr!i2)2 + (2ftr!23)2 + A2, 



(4) 



Our aim is to induce population transfer from state 
1) to |3) by maintaining the system in state |Po} and 
changing the characteristics of jPo) from |1) at t = to 
3) at t — ^max by appropriate control of the tunnelling 
rates, without population leakage into the other eigen- 
states. Provided we start in an eigenstate, and transform 
our Hamiltonian adiabatically, we will remain in the same 
eigenstate. The criterion for adiabaticity is jS| 



1^0 -^±1 » l(I>o|2?±)|- 



(5) 



It is usual to apply Gaussian pulses to realize the adi- 
abatic transfer, and we consider pulses of the form 



a(3 



(6) 



FIG. 2; Energies of the states |E>+), |I?o> and \V-) from top 
to bottom, respectively, through the CTAP pulse sequence 
described in Eq. [7|with A = 0. 



where i^/S and Uap are the peak time and standard devi- 
ation of the control pulse modulating the tunnelling rate 
between states \a) and \[3). For simplicity we set the max- 
imum tunnelling rates and standard deviations for each 
transition to be equal, i.e. il™^'^ = and aap = cr, 

and set A = in the discussion that follows. Under 
these conditions, transfer is optimized when the width of 
the pulses equals the time delay between the pulses [l9| . 
with total pulse time imax, we choose a = tmax/S and the 
barrier controlled tunelling rates 



ni2{t) = l^'^'^'^exp 
VL2z{t) = f^'^^'^exp 



(7) 



This ordering, where f223 is applied before f2i2 is known 
as the counter-intuitive pulse sequence and has significant 
advantages in improving transfer fidelity over other pulse 
sequences. Under these conditions, the energies of the 
eigenstates are shown in Fig.[51for tmax = lOir/Vl. 

In order to proceed, we numerically solve the master 
equations for the density matrix, p. 



r [p - diag(p)] , 



(8) 



where T is the T2 (pure dephasing) rate, assumed to act 
equally on all coherences. As we are primarily consid- 
ering solid-state systems, we ignore the Ti rates, as T2 
is generally much faster, although such rates could be 
added in principle. 

In Fig. 121 we present results showing the population 
transfer using the counter-intuitive and intuitive pulse 
orderings without dephasing. As long as the adiabaticity 
criteria is satisfied, the fidelity of population transfer is 
very high. In practise the maximum possible transfer 
rates will be a few times greater than Vl^'^^. 

We have considered the effects of T2 dephasing in Fig.0| 
which shows the populations as a function of the to- 
tal transfer time, imax and F. We can broadly identify 
three regimes here. When t„i!aS^™^^ / S 2, nonadiabatic 
transfer is observed, typified by the oscillating popula- 
tions in jl) and |2). Note that no appreciable population 
builds up in |3) in this regime. For tmnK^'^'^^ / tt > 5 and 
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FIG. 3: Charge transfer using both counter-intuitive 
[Fig-EIa) -OIc)] and intuitive [Fig.|3d) -Elf)] pulse sequences 
with no dephasing (F = 0). (a) Tunnelling rates as a func- 
tion of time (in units of vr/Sl™^", S7i2 is the solid line, Q23 is 
the dashed line, (b) Populations as a function of time as the 
trajectory of (a) is followed. Note the smooth migration of 
population from |1) to |3) with no population in |2). (c) Final 
populations as a function of tmax (in units of tt/SI™^'') where 
the pulse sequence and pulse widths are determined by Eq.|7| 
(d) Intuitive pulse sequence, as before Q12 is the solid line and 
f223 is the dashed line, (e) Populations as a function of time 
along the trajectory of (d). In this case there are significant 
oscillations and a relatively large population in |2), this pulse 
ordering is clearly unsatisfactory for high fidelity transport, 
(f) as (c) but for the intuitive pulse ordering. Note that in 
this case, to get high-fidelity transfer, more time is required. 



r/il"^^^ > O.I the populations of all three states tend to 
equilibrate, again this regime is not useful for high fidelity 
transfer. In order to maintain high fidelity, adiabatic 
transfer, we need tm^M'^'''' /n « 3.75 and F/fJ™^'^ < O.I. 

That there is a window of transfer times to achieve 
high fidelity transfer should not be surprising, consid- 
ering that the inclusion of T2 effects introduces a char- 
acteristic timescale, beyond which coherence cannot be 
maintained. So one must perform a transfer slow enough 
to be adiabatic, but fast enough to not be decohered. Al- 
though these values for T/il'^^^ may appear small com- 
pared with optical systems, one should realise that in 
STIRAP we are comparing the transfer times with the 
two-photon decoherence of the two ground states, which 
is normally very much longer than the decoherence rates 
of the optical transitions. In our case wc have assumed 
that the action of the environment is to decohere all su- 
perposition states equally as a pure dephasing process. 
Pure dephasing is known to limit the performance of 
charge qubits and is manifested by electrostatic inter- 
actions with background charge fluctuators and phonon 
interactions 9, 20]. 

Because CTAP is a coherent transfer method, there 
will be no heat dissipated by the electron during the 
transfer, and because the initial and final states are 
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FIG. 4: Pseudo-color plots showing populations (scale: red — 
1, blue = 0) (a) pn, (b) p22 and (c) P33 as a function of total 
pulse time (in units of Tr/^Y^^^^) and dephasing, F (in units of 
j-jmax-j Successful transfer is identified by pn ~ P22 ~ and 
Pas ~ 1, this condition is satisfied when F/fi'"^^ 0.1, and 

tmaxJ^™""/^ ~ 3.75. 




FIG. 5: Pseudo-color plot with overlayed contours showing 
P33 as a function of energy mismatches A31 and A21 in units 
of n™"" for tniax = 157r/r2™'"'. The system is less sensitive to 
A21 than A31 as expected for three-state adiabatic passage. 



experimental parameters are known to arbitrary preci- 
sion. In general, however, the characterization of exper- 
imental Hamiltonians is a difhcult problem, and com- 
plicated protocols must often be employed in order to 
extract the parameters to high fidelity (see for exam- 
ple Schirmer et al. |0| for the case of a two-state sys- 
tem). One advantage of STIRAP like schemes is that 
they are relatively insensitive to the exact conditions. 
We have illustrated this in Fig. El where the effects of en- 
ergy mismatches between the states have been modeled. 
Here we show P33 as a function of A31 = E3 — Ei and 
A21 ^ E2 — E1 for tinax — 157r/f2 to simulate the effect of 
a systematic error in the energy levels. We have ignored 
the effects of dephasing for clarity in this figure. As is 
expected from the form of the dressed states in Eq.[21the 
behaviour is very insensitive to A21 but more strongly 
influenced by A31. Note that an energy mismatch of as 
much as r2™'*'^/5 still permits P33 ~ 0.9. One could imag- 
ine optimizing the transfer conditions by studying the P33 
and feeding results back to 'tweak' the energies. 



TRANSPORT THROUGH A MULTI-DOT 
SYSTEM 



One obvious extension of this work is to consider the 
passage through more than one intervening dot. This 
extension has been considered several times in the optical 
literature and we concentrate here on one particular 
extension, the so-called straddling scheme of Malinovsky 
and Tannor |22l |. For simplicity we will not include the 
effects of dephasing in this treatment. 

To realise the straddling tunnelling sequence (SCTAP) 
we must augment the original pulse sequence of Eq.|Z|by 
the straddling pulses which are the same for all interven- 
ing tunnelling rates. To be more precise, if we label the 
dots 1,2 •■•n with tunnelling rates r2i2, $723 ■ • ■ i^n-i,n 
and we wish to transfer population from dot 1 to dot 
n, then we apply the pulse sequence 



r!i2(i) = r!i = r!grexp<^- 



[t- (W/2 + g)]' 
2^2 



ground states, there can be no heat dissipation after 
transfer either. Therefore, provided the transfer is fast 
relative to the decoherence time, there will be no heat 
dissipation regardless of the transfer speed. This should 
be contrasted with the quasi-adiabatic result for QDCA 
switching where there is a necessary tradeoff between 
switching time and heat dissipation [l^ . It must be 
stressed that although our scheme does not dissipate en- 
ergy in the electronic system, there will of course be en- 
ergy dissipation involved with charging/discharging the 
control gates. 

In the preceding discussion, we have assumed that the 



n„-i,„W = f72 = f7crexp' 



[t - (imax/2 - CJ)Y 



r!,,,+i(t) = f7s = f^r'exp 



2a2 

{t tmax/2) 

4^^ 



,(9) 



where 1 < z < 71— 1, VL™f^ is the maximum tunnelling rate 
for the counter-intuitive pulses, and Vl™^^ is the maxi- 
mum tunnelling rate for the 'straddled' transitions. Note 
that in contrast to the 3 state pulse sequence in Eq. 
(with pulse separation a), the spacing between pulses is 
2ct for the straddling scheme. The pulse separation in 
each case was chosen to optimize the transfer fidelity. 
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FIG. 6: Energies of the 5 (a) and 6 (b) dot systems with 
pulse sequences defined by Eq. |^and tmax = lOyr/n™/'' and 
^max ^ sf^gf. Notice that the 6-state system has, in fact, 
two central states (with neither at zero energy) whereas the 
5-state system has one central state at zero energy. The two 
central states of the 6-state system are highlighted in the in- 
set, which is an expanded view of the boxed region about the 
centre (the horizontal scale has been omitted for clarity, and 
is [3-7]^/ngr)- 



We are unable to provide a simple physical explanation 
for why the pulses separations in each case are different. 

The Hamiltonian for the straddling scheme is (A = 0) 
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f7s 
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172 



(10) 



In contrast to Malinovsky and Tannor [2j|, we find that 
this scheme is only effective for odd numbers of dots. 
This discrepancy between even and odd numbers of states 
can be understood by examining the eigenvalue structure 
through an adiabatic passage. 

In Fig. we show the eigenvalues for the 5 and 6 
state SCTAP pulses with = bn'^f and t^ax = 

107r/f2™j'^. Note that for 5 state transfer, there is one 
state at zero energy, with the other states symmetric 
about it, whereas for the 6 state system, there are two 
states close to zero. These patterns are typical of all 
even and odd state systems. The proximity of these two 
states to zero means that the adiabaticity criterion is 
much harder to satisfy for states with an even, rather 
than odd number of states. 

The eigenvalues of the Hamiltonian in Eq. 1101 are in 
general rather complicated, however, the zero energy 
eigenstate has a relatively simple form in this case which 



roo 

0.99 
0.98 
0.97 
0.96 
\ 0.95 
0.94 
0.93 
0.92 
0.91 
0.90 







5 


/ 11 


r 


/ 

































22 



28 



FIG. 7: pnn as a function of fmax (in units of yr/f^™/'') for (l^ft 
to right along tmax axis) 5,7,9,11,13, and 15 state systems for 
constant fig""' = imcT- 



is (unnormalized) 

^ cosOill) - (-l)"sinei|2n- 1) 



X 



n-l 

i=2 



(11) 



where Qi is defined as above, and 



X = 



(12) 



for an array of 2n — 1 dots. We can now see that the 
scheme prevents occupation of the even numbered states, 
and minimizes population in the other undesirable states, 
as long a,s X ^ 1. Note the similarities between these 
states, and the one s derived in the context of the chain A 
scheme in Ref. ^231 which had alternating couplings and 
energy detunings. 

To determine the parameter range needed to achieve 
high fidelity transfer, we solve the density matrix equa- 
tions of motion for the multi-dot system, with varying 
tmax and for constant f^g^^'^ = 50. In Fig. H we 

present results showing p„„ as a function of the total 
transfer time, imax for 5,7,9,11,13, and 15 state systems. 
Note the almost monotonic behaviour with occasional 
'ripples'. Because of the different pulse separations of the 
3 state case, it is not valid to directly compare it with the 
other cases, although population transfer for the 3 state 
case is always easier to achieve than for other cases. 



QUANTUM INFORMATION TRANSFER 
ACROSS A LINEAR CHAIN 



We now consider applying the CTAP scheme to the 
problem of transferring quantum information across the 
chain. This may be of use in quantum computation, 
where quantum information is required at a different 
point from where it is generated (e.g. taking informa- 
tion from a 'computation zone' to a 'measurement zone'). 
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It is important to realise that the CTAP schemes as de- 
scribed are not capable of realizing swap operations with- 
out some modification, perhaps along the lines suggested 
by Unanyan et al. |24l| . However it is possible to perform 
a transfer of quantum information across a network. We 
consider the problem of attempting to switch the super- 
position state of a single charge-qubit, where the qubit 
is defined by a single electron in a superposition of two 
sites, IV'} = +/?|2} with unoccupied sites |3) and |4). 
We wish to transfer the complete quantum state in |2) to 
|4) leaving the rest of qubit unchanged. Such information 
transfer would correspond to the following transition: 



case IS 



a|l)-H/3|2)^a|l)+/3|4). 



(13) 



One can show quite easily that such a transformation can 
be achieved by simply applying the standard CTAP pulse 
sequence to the 2, 3, 4 system, without affecting state 
The Hamiltonian of the operation is 
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(14) 



7i = 



and the initial and final states are 
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(17) 



l^,)=a|ii)+/3|i?i), 

|V/)=«|i3)+/3|i?3), 



(18) 



respectively. Note that in this case the tt rotations ob- 
served earlier have effectively cancelled so that the final 
state is in phase with the initial state. 

One may extend these calculations to transfer across 
many dots using an extension of the SCTAP sequence. 
However for quantum information transfer, transfer 
across an even number of dots will be unsatisfactory as all 
states change their energy along the SCTAP trajectory. 



CONCLUSIONS 



with fli and fl2 defined by fli2 and 0,23 from Eq. Irre- 
spectively. If we start with our system in an arbitrary 
superposition state 

= a\l) + P\2) (15) 

Numerically solving the master equations along the 
CTAP pulse trajectories yields the final state 

IV'/) =a|l)+e*"/3|4) (16) 

which is the desired final state up to a constant tt phase 
rotation, which can be understood by inspection of the 
eigenstates in Eq. |21 Wc note that without maintaining 
the zero energy condition of the passage state (in this 
case the equivalent of the |Po) state from Eq.^J, the sys- 
tem would pick up an extra phase term due to the energy 
difference between state |1) and the adiabatically trans- 
formed state which depends critically on the exact tra- 
jectory taken. Such a consideration clearly shows the im- 
portance of adiabatic transform mechanisms where there 
is no change in energy and are therefore more robust 
against errors and fluctuations. 

We may also consider the case where we wish to trans- 
fer both parts of the qubit using CTAP pulses. Consider 
first the chain |L3), |L2), |L1), |i?l), |i?2), |i?3} where 
we desire to transfer an arbitrary superposition of \L1) 
and \R1) to \L3) and |i?3). In this case we can apply the 
CTAP pulses simultaneously. The Hamiltonian in this 



We have described a method of coherent electronic 
transport through a triple-well system which we term 
CTAP (Coherent Tunnelling Adiabatic Passage) by anal- 
ogy with the optical process STIRAP (Stimulated Raman 
Adiabatic Passage) . This scheme is realized by electronic 
modulation coherent tunnelling rates by control of bar- 
rier gates. It is a high fidelity process, provided the co- 
herent tunnelling rates are fast compared with the T2 
dephasing times. This scheme involves no change in the 
electronic energy throughout the transfer regardless of 
transfer speed, and so provides an interesting alternative 
to more conventional methods for quasi-adiabatic clock- 
ing in QDCA cells. 
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